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p^ ABSTRACT 

(-H This paper presents an environment and stellar mass study of a large sample of star- 

O . forming Ha emitters at z = 0.84 from the High-z Emission Line Survey (HiZELS), over 

I 1.3 deg2 spHt over two fields (COSMOS and UKIDSS UDS). By taking advantage of a ti-uly 

2 panoramic coverage of a wide range of environments, from the field to a rich cluster, it is 

'(-J shown that both stellar mass and environment play crucial roles in determining the properties 

^ of star-forming galaxies. Specific star formation rates (sSFRs) decline with stellar mass in all 

i__i environments, and the fraction of Ha star-forming galaxies declines sharply from « 40 per 

cent for galaxies with masses around 10^° M© to effectively zero above 10^^'^ M©, confirm- 

Cn ing that mass-downsizing is generally in place by z ^ 1. The fraction of star-forming galaxies 

>^ is also found to fall sharply as a function of local environmental density from w 40 per cent in 

^_S the field to approaching zero at rich group/cluster densities. When star formation does occur 

^^ in such high density regions, it is found to be mostly dominated by potential mergers and, 

indeed, if only non-merging star-forming galaxies are considered, then the environment and 

mass trends are even stronger and are qualitatively similar at all masses and environments, re- 

r^ spectively, as in the local Universe. The median star-formation rate of Ha emitters at z = 0.84 

is found to increase with density for both field and intermediate (group or cluster outskirts) 
densities; this is clearly seen as a change in the faint-end slope of the Ha luminosity func- 
tion from steep (a w —1.9), in poor fields, to shallow (a sa —1.1) in groups and clusters. 
Interestingly, the relation between median SFR and environment is only found for low to 
moderate-mass galaxies (with stellar masses below about lO^^'^M©), and is not seen for the 
most massive star-forming galaxies. Overall, these observations provide a detailed view over 
5^ a sufficiently large range of mass and environment to reconcile previous observational claims: 

stellar mass is the primary predictor of star-formation activity at z ^ 1, but the environment, 
while initially enhancing the median SFR of (lower-mass) star-forming galaxies, is ultimately 
responsible for suppressing star-formation activity in all galaxies above surface densities of 
10-30 Mpc^^ (group and cluster environments). 

Key words: galaxies: high-redshift, galaxies: luminosity function, cosmology: observations, 
galaxies: evolution. 



1 INTRODUCTION 

At the present epoch, star formation activity is strongly depen- 
dent on environment (e.g. Lewis et al.||2002[ [Gomez et al._2003{ 
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Based on observations obtained using the Wide Field CAMera on the 



[Tanaka et al.|2004{ [Mahajan et al.|20To i: while clusters of galax 



3.8m United Kingdom Infrared Telescope, as part of the High-^ Emission i^s seem to be primarily populated by passively-evolving galax- 

Line Survey (HiZELS) i^S' star-forming galaxies are mainly found in low-density envi- 

t E-mail: drss@roe.ac.uk ronments l[Dressler|1980^. It is well-established (e.g.|G6mez et al.| 
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[2003i|Kauffmann et al.|2004)|Best|20"04l l that the typical star for- 
mation rates of galaxies - and the star-forming fraction - decrease 
with local galaxy density (often projected local density, E) both 
in the local Universe and at moderate redshift (e.g. Kodama et al.| 
[2004l l. Indeed, at « ~ 0.2 - 0.3, studies such as |Couchet al.H2001> 
or |Balogh et al.|f2002| l found that the Ha luminosity function in 
rich clusters, whilst having roughly the same shape as in the field, 
have a much lower normalisation (~ 50 per cent), consistent with 
a significant suppression of star formation in such environments. 

Typically, active star-forming galaxies in the local Universe 
are also found to have lower masses than passive galaxies and, in- 
deed, the most massive galaxies are mostly non-star-forming (an 
observational result often known as mass-downsizing; |Cowie et al.| 
|1996j ); this implies that the most massive galaxies assembled their 
stellar mass quicker and have not formed a significant amount of 
stars since then. Massive galaxies are predominantly found in high 
density environments, but the mass-downsizing trend is not simply 
a consequence of the enviroiraiental dependence, nor vice-versa. 
Recent studies such as [Mahajan et al.| ^2010^ - and references 
therein - show that the environment trends hold for both massive 
and dwarf galaxies, and that the mass trends also hold for both field 
and cluster environments. [Peng et al.|p010j also finds that the ef- 
fects of mass and environment are completely separable. 

When did the environment and mass dependences of star- 
forming galaxies start to be observable, and how did they affect 
the evolution of galaxies and clusters? Clearly, in order to prop- 
erly answer such questions it is mandatory to conduct observational 
surveys at high redshift, which can then be used to test theoretical 
models of galaxy evolution. By z ~ 1, some authors have claimed 
to have found a flattening or even a definitive reverse of the rela- 
tion between star formation activity and local galaxy density, either 
by studying how the average/median star formation rates of galax- 
ies change with local density (e.g. plbaz et al.|2007[ [Cooper et al.] 
|2008 ) or by looking at the star-forming fraction as a function of 
density (e.g. Ideue et al.|2009f. Similar trends have also been sug- 
gested at 2 ~ 2 by Kodama et al.| ^2007^ , and recently |Hayashi| 
|et al.|pOTo) found that the fraction of star-forming galaxies seems 
to remain constant up to the highest densities in a cluster environ- 
ment at z ~ 1.5. These results would be naturally interpreted as 
a sign of a clear evolution if other studies (e.g. Finn et al._2005, 
[Patel et al.||2009[ [Koyama et"Sr||2010| ) had not indicated the op- 
posite result. The latter studies argue that even at « ~ 1, and for 
the richest/cluster environments, both star formation rate and the 
star-forming fraction decline with increasing local density. Never- 
theless, Koyama et al.| ( |2010^ also found that the Hq luminosity 
function of a rich cluster at z ~ 0.8 has roughly the same shape 
as the field luminosity function from |Sobral et al.H2009a[ >, but with 
a higher normalisation; this implies that even if the environment is 
already suppressing star formation by z ~ 1, there is still a consid- 
erable population of star-forming galaxies in high-redshift clusters 
(c.f. [Butcher & OemIer| 19841 1. 

Part of the discrepancies in environmental trends between dif- 
ferent studies may well be connected with possible mass depen- 
dences already in place at high redshift. While at z ~ the large 
and deep surveys can now probe very different mass regimes with 
large enough samples within each environment to confirm and dis- 
tinguish mass and environment trends individually (c.f. |Peng et al.[ 
[2010j ), that is clearly not the case at high redshift, where luminosity 
limits inevitably bias samples towards higher mass galaxies. Thus, 
if mass-downsizing is already in place at high redshift, and mas- 
sive galaxies are preferentially found in high density regions, sam- 
ples with different luminosity limits and/or based on different se- 



lections may find contradictory results. Indeed, such downsizing 
trends have now been identified out to z ~ 1.5, by revealing that 
the most massive galaxies are undergoing very little or no star for- 
mation, or by finding a population of massive galaxies already in 

[Marchesini et al. 



the red sequence (e.g. Bauer et al. 



[Santini et al.[[2009l [Fontanot et al.||2009| [Cirasuolo et al.||2010l 
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Therefore, in order to identify and distinguish the separate roles of 
mass and environment on star formation at high redshift, one really 
requires clean, robust and large samples of star-forming galaxies 
residing in a wide range of envirormients and with a wide range of 
masses, together with samples of the underlying population found 
at the same redshift. 



Narrow-band Ha surveys are one of the most effective ways 
to gather representative samples of star-forming galaxies at differ- 
ent epochs, and the scientific potential of these is now being widely 
explored, following the development of wide-field cameras in the 
near-infrared. This has led to significant progress (e.g. Ly et al.| 



2007: Shioya et al. 2008; Geach et al. 'loSSJ [Sobral et al.[[2009af 
Garn et al. 2010; Sobral et al. 2010) taking advantage of large 



samples of typically hundreds of Ha emitters up to z = 2.23. 
In particular, HiZELS, the Hi-Redshift(z) Emission Line Survey 
( [Geach eFa l. 2008, S obral et al. [[2009a} , is playing a fundamen- 
tal role by obtaining and exploring unique samples of star-forming 
galaxies selected through their Ha emission (and other major emis- 
sion lines; [Sobral et al.[[2009bl l, redshifted into the J, H and K 
bands (see Best et ar][2010 1 for an overview). By using a set of 
existing and custom-made narrow-band filters, HiZELS is survey- 
ing Ha emitters at z = 0.84, z = 1.47 and z = 2.23 over several 
square degree areas of extragalactic sky. Narrow-band surveys such 
as HiZELS probe remarkably thin redshift slices (Az ~ 0.03), 
allowing the study of very well-defined cosmic epochs. They sur- 
vey the redshift slices down to a known flux limit and they can se- 
lect equivalent populations at different redshifts. Combined with a 
wealth of high-quality multi-wavelength data in the targeted fields, 
wide narrow-band surveys like HiZELS can be fully explored to 
detail the roles of "nature" and "nurture" in galaxy evolution, by 
unveiling the importance of the environment on a well-defined star- 
forming population at a single epoch. 



This paper presents a detailed study of the dependencies of 
star formation on environment and mass, and their relations with 
galaxy morphology, using a large sample of narrow-band selected 
Ha emitters at z = 0.84. Using this sample, it is possible to recon- 
cile many of the apparently contradictory results in the literature. 
The paper is organised in the following way. §2 outlines the data 
and selection, the samples, their sky distribution and the mass and 
local density estimations. §3 presents the dependence of star forma- 
tion activity on mass, together with the downsizing trends found at 
z = 0.84. §4 details the connection between star formation and 
the environment, presenting the relations between local density, 
star formation activity and star-forming fraction for Ha emitters 
at z = 0.84. §5 presents a detailed, combined mass-environment 
view of star formation at z ~ 1. §6 discusses and interprets the re- 
sults in the context of other studies and finally §7 gives a summary 
and the conclusions. An Ho = 70kms^^ Mpc^^, Q,m = 0.3 and 
f^A = 0.7 cosmology is used and, except where otherwise noted, 
magnitudes are presented in the AB system. 
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2 SAMPLES AND PROPERTIES 

2.1 The sample of robust Ha emitters at z — 0.84 

This study uses the large sample of Ha emitters at z = 0.84 from 
HiZELS presented in Sobral et al. (2009a, S09a). Briefly, the sam- 
ple was derived from a narrow-band J filter (Acfi = 1.211(j.m, 
AA = 0.015^m) survey using the Wide Field CAMera (WF- 
CAM) on the United Kingdom Infrared Telescope (UKIRT), reach- 
ing a flux limit of 8 x 10~^^ergs~^ cm~^ over ~ 0.7 deg^ in the 
UKIDSS Ultra Deep Survey (UDS, [Lawrence et"ar][2007l l field, 
and ~ 0.8 deg^ in the Cosmological Evolution Survey (COSMOS, 
[Scoville et al.|2007| field. The survey resulted in the selection of 
1517 line emitters which were clearly detected in NBj (signal- 
to-noise ratio > 3) with a J-NBj colour excess significance of 
E > 2.5 and observed equivalent width EW> 50 A - see S09a for 
full details on the selection. 

Photometric redshiftaj were used to select a sample of 743 
Ha emitters over a co-moving volume of 1.8 x 10^ Mpc'^ at z — 
0.845. Of these, 477 are in COSMOS (0.76 deg^) and 266 in UDS 
(0.54 deg^). The completeness and reliability of this sample were 
studied using the ~ 10* available redshifts from 2;C0SM0S Data 
Release 2 jLilly et al.|2009, ), spectroscopically confirming ~ 100 
Ha emitters within the photometric redshift selected sample; this 
allowed to estimate a > 95 per cent reliability and > 96 per cent 
completeness for the sample in COSMOS. Here, the sample of Ha 
emitters is modified from that in S09a in three ways: 
(i) Following [Sobral et al.|J20r0j , the Ha sample is slightly mod- 
ified on the basis of the new photometric redshifts, incorporating 
a higher number of bands, deeper data and accounting for possi- 
ble emission-line flux contamination of the broad-bands (c.f. [Ilbertl 
[et al.[2009l Cirasuolo et al. in prep.). In particular, the revised sam- 
ple in UDS contains 257 Ha emitters over 0.52 deg^, whilst that in 
COSMOS is unchanged (477 emitters over 0.76 deg^). 
(ii) Some of the Ha emitters may contain an AGN, which can be 
responsible for a significant fraction of the Ha flux and are capable 
of changing the host galaxy spectral energy distribution. Within the 
Ha sample, AGN candidates have been identified by [Gam et al.[ 
( [20 10^ using a wide range of methods. The analysis led to the iden- 
tification of a maximum of 74 AGN (40 in COSMOS and 34 in 
UDS). From these, 40 are classified as possible AGN and the re- 
maining 34 as likely AGN - implying 5—11 per cent AGN contam- 
ination of the Ha sample at z — 0.84. For the analysis presented in 
this paper all 74 AGN candidates (possible and likely) are excluded 
from the star-forming Ha sample (though they are kept in the un- 
derlying population sample), but it has been checked that if these 
are retained then none of the results or conclusions of this paper 
are altered. This results in a robust sample of 660 star-forming Ha 
emitters at 2; — 0.84. 

(ill) Finally, as detailed in §2.3, the sample is also modified to in- 
clude potential Ha emitters with low equivalent widths. 



available. These include deep intermediate and broad-band imag- 
ing from the ultra-violet (UV) to the infrared (IR), making it possi- 
ble for excellent photometric redshifts (photo-zs) to be determined. 
These data are fundamental to test downsizing and environmental 
trends, since these trends can only be investigated by selecting an 
underlying population of galaxies within the same redshift range as 
the Ha population. Whilst the underlying population will never be 
as clean as the Ha sample, the high quality photometric redshifts at 
z ~ 0.8 available in COSMOS and UDS are good enough to do a 
first-pass cut at emulating the narrow-band filter selection. For the 
COSMOS field, the zCOSMOS DR2 secure (> 99 per cent confi- 
dence) spectroscopic redshifts can then be used to study the com- 
pleteness and the contamination of various photometric-redshift se- 
lected samples (as well as to improve both, by taking out confirmed 
contaminants and introducing spectroscopically confirmed galax- 
ies missed by the photometric redshifts). These results provide ap- 
propriate corrections to the photometric redshift selection. In Ap- 
pendix A, a range of different initial photo-z cuts are considered 
and the resultant correction factors based on the contamination and 
completeness are calculated; these are used to show that the results 
in the paper are robust to changes in the photometric redshift cut. 

The analysis in this paper is carried out using the background 
population sample with 0.82 < Zphoto < 0.87 for COSMOS 
(where photo-zs are very precise; the sample is also corrected and 
improved using spectroscopic redshifts), estimated to be 77 per cent 
complete and with a contamination of 59 per cent (when compared 
with the redshift range covered by the narrow-band top hat filter 
profile, 0.83 < z < 0.86). For UDS, the underlying sample is 
obtained with galaxies having 0.79 < Zphoto < 0.9. Reliable con- 
tamination and completeness estimates are not available due to the 
lack of large spectroscopic data-sets, but the completeness can still 
be roughly estimated based on the recovery of Ha emitters (see 
Appendix A), yielding ~ 60 per cent. Also, given the uncertainty 
on the photo-2S, a ~ 75 per cent contamination is expected and 
this is assumed throughout the analysis. This, of course, results in 
a higher uncertainty in the following analysis for UDS, but as the 
next sections will fully detail, similar trends are recovered both in 
COSMOS and UDS. 

A common K < 23 (observed, corresponding to rest-frame 
J at z = 0.84) cut is applied (prior to the completeness and con- 
tamination studies) to match the K completeness of the Ha emit- 
ters sample and minimize potential biases towards lower luminos- 
ity, lower stellar mass galaxies; this also guarantees that both the 
Cirasuolo et al. (in prep.) and the [Ilbert et al.[p00 9') catalogues are 
more than 90 per cent complete and that the photometric redshifts 
perform to their best at the studied redshift, avoiding introducing 
extra biases into the analysis. Sources clearly classified as stars are 
also rejected. The final underlying sample in COSMOS contains 
3656 sources (~ 4800 per square degree) including all Ha emit- 
ters, while the equivalent sample in UDS contains 2688 galaxies 
(~ 5100 per square degree). 



2.2 The underlying population samples 

As the samples were obtained in two of the best-studied square de- 
gree areas of equatorial sky, a wealth of multi-wavelength data are 



^ Photometric redshifts used in S09a for COSMOS present a(Az) = 
0.03, where Az = (zpi^ot — 2spec)/(l + ^spec)', the fraction of outliers, 
defined as sources with Az > 3cr(A2), is lower than 3 per cent, while for 
UDS, the photometric redshifts have a{Az) = 0.04, with 2 per cent of 
outliers. 



2.3 The potential Ha emitters at z = 0.84 with EW< 50 A 

The equivalent width (EW) limit was adopted by S09a to improve 
the selection of emission-line objects and, in particular. Ha emit- 
ters at z — 0.84. It helps to avoid contamination both from bright 
sources, which can have an apparently significant narrow-band ex- 
cess due to small colour terms, and from weak (low EW) emission 
lines, which may be close in wavelength to Ha and thus difficult to 
distinguish with photo-2S. However, the EW limit can lead to ne- 
glecting galaxies which have emission-line fluxes above the limit of 
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the survey but have a relatively luminous stellar continuum in the 
J band (rest-frame Rai z = 0.84). In S09a, this is dealt with by 
conducting detailed simulations which replicate the selection crite- 
ria. Simulations are used to compute a completeness correction as 
a function of luminosity for deriving the Ha luminosity function 
(S09a). These estimate that a total of « 10 per cent of all Ha emit- 
ters above the flux limit may be missed by the EW> 50 A selection 
in both COSMOS and UDS. However, they do not provide infor- 
mation about how this varies as a function of mass or environment. 

Motivated by the need to provide a fair comparison of the Ha 
selected sample with the underlying population at the same red- 
shift, the latter samples are used (in COSMOS and UDS). For 
all photo-z selected galaxies, data from S09a are used and the 
EW > 50 A requirement is waived. Galaxies within the photo- 2; 
sample that present emission line fluxes (colour excesses S > 2.5, 
following S09a) above the survey limit are considered to be poten- 
tial additional low EW Ha emitters. 

This approach results in selecting 76 additional sources 
in COSMOS and 42 in UDS above a flux hmit of 8 x 
10~^'^ergs~^ cm~^ (these sources have EW= 35 ± 11 A, with 
the lowest EW being 13 A). For COSMOS, 25 of the 76 recov- 
ered sources have high quality spectroscopic redshifts and 17 are 
indeed confirmed to be Ha emitters at z — Q.84~| implying a con- 
tamination of 32 per cent (clearly demonstrating the initial motiva- 
tion for the EW limit). After removing the confirmed contaminants, 
the sample is reduced to 68 Ha emitters. Furthermore, inspection 
of ^COSMOS spectra of the confirmed Ha emitters reveals [Oil] 
emission lines with fluxes ~ 1 x lO^^^ergs"^ cm~^, consistent 
with the estimated weak Ha emission for these sources. None of the 
sources present H/3 or [Olll] with S/N > 1, although 2 of these 
sources present a red, detectable continuum, with strong absorp- 
tion lines such as H and K. Only one source is detected as a strong 
24 |i.m source; this is flagged as an AGN in the general COSMOS 
catalogue. 

For the UDS field, the lack of spectroscopic data does not al- 
low for a similar investigation to be conducted, and with a wider 
photo-z cut being used, a 40 per cent contamination (by both lines 
close enough to Ha and others, as found in COSMOS) is assumed 
for the low EW Ha sample in UDS. 

A sample containing a total of 770 Ha emitters will be used in 
this paper, but the estimated contamination rates of the samples of 
potential Ha emitters are taken into account by using those as rel- 
ative weights in the analysis. Ha luminosities and SFRs for all Ha 
emitters are derived following S09a, accounting for dust extinction 
as in Gam et al. (2010, equation 10), but it should be noted that the 
results remain unchanged if no, or a uniform extinction is applied 
to the sample. 



2.4 Morphological classification with HST imaging 

The COSMOS field has sensitive Hubble Space Telescope (HST) 
ACS F814W coverage, providing detailed morphological informa- 
tion for all samples. Automated classifications obtained by ZEST 
(Scarlata et al. 2007) are pubhcly available for the entire COSMOS 
field, which classifies galaxies in 3 morphological classes: 1) Early- 
types, 2) Disks/spirals and 3) Irregulars. Reliable estimates exist 
for 88 per cent of the complete underlying sample in COSMOS at 



^ The data also reveal 1 [Arm] emitter aX z = 0.69, 2 [SlI] emitters at 
z = 0.79, 1 [Nil] emitter (z = 0.82), and 3 [Ol] emitters at 2 = 0.89, 
while 1 source is not identified (z Ri 0.99). 



Table 1. Parameters used to generate the SED templates using the Bruzual 
& Chariot (2003) package and the new Bruzual (2007) models. 
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0.82 < Zy < 0.87. Visual classifications were also obtained by 
S09a for Ha emitters at z = 0.84; the authors find that the two 
agree very well, with only small differences arising from the use 
of colour information in S09a (from Subaru data), while ZEST uses 
the HST F814W imaging only. 

Due to a considerable number of sources showing evidence of 
potential merging activity in S09a, the Ha sample was also clas- 
sified independently into merger classes. Data from Subaru (B, r' 
and z) was combined with the resolution from HST so that the 
colour information could be used to obtain a good identification of 
potential mergers (e.g. distinguishing between normal disks with 
strong dust bands and real mergers). A source was classified as a 
potential merger when: 1) it presented a clearly disturbed morphol- 
ogy or disturbed disk, inconsistent with being that of a normal disk 
galaxy; or 2) it presented more than one bright point-like source and 
the colour information was inconsistent with one of those being a 
spiral arm; or 3) there were two or more galaxies which were very 
close (< 15kpc) - see S09a for examples of each morphological 
class. Here, the same criteria are used to extend this classification 
to galaxies in the sample of potential Ha emitters with EW< 50 A. 
These sources are found to be mostly discs (65 per cent), with 15 
per cent irregulars and 20 per cent early-type galaxies (the signifi- 
cant increase in the fraction of early-type galaxies is the most no- 
ticeable difference between the populations with different EWs, as 
the sample with lower EWs is able to probe much lower specific 
star formation rates); 33 per cent of the sample shows evidence of 
potential merging activity. 

2.5 Stellar mass estimates 

Deep multi-wavelength data are available for all Ha emitters and 
galaxies within the underlying populations and these data are used 
to compute stellar masses. The multi-wavelength data are used to 
perform a full SED x^ fit - normalised to one solar mass - to each 
galaxy; the stellar-mass is the factor needed to re-scale the lumi- 
nosities in all bands from the best model to match the observed 
data. This method is similar to that used by [Fontana et al.| (2004b 
and |Ilbert et al.|(2010^ . The SED templates are generated with the 
stellar population synthesis package developed by 'Bruzual & Char^ 
lot (2003 ), but the models are drawn from Bmzual ( 2007 1. SEDs 
are produced assuming a universal initial mass function (IMF) from 
|Chabrier] ( [2003[ l and an exponentially declining star formation his- 
tory with the form e~'' ^, with r in the range 0.1 Gyrs to 10 Gyrs. 
The SEDs were generated for a logarithmic grid of 200 ages (from 
0.1 Myr to 6.5 Gyr ~ the maximum age at z = 0.84). Dust extinc- 
tion was applied to the templates using the |Calzetti et al.|(2000[ l law 
with E{B — V) in the range to 0.5 (in steps of 0.1). The models 
are generated with three different metallicities, including solar; the 
reader is referred to Tablefllfor a list of the parameters used. 
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For the SED fitting, it is assumed that all Ha emitters are at 
z = 0.845 (the mean redshift for detecting Ha with the narrow- 
band filterjjand the full filter profiles are convolved with the gen- 
erated SEDs for a direct comparison with the observed total fluxes. 
Up to 36 wide, medium and narrow bands (from the GALEX far 
ultra-violet and near ultra-violet bands to Spitzer 4 IRAC bands, c.f. 
[Ilbert et al.|2009l|20T0) are used for the COSMOS field, after ap- 
plying all of the necessary corrections to obtain the total fluxes, as 
suggested and fully detailed in Capak et al. (2007) and |Ilbert et al.| 
( |2009[[2010[ l: individual aperture corrections, average systematic 
offsets per band and individual galactic extinction per band. For 
the UDS field, fewer bands are available (medium band imaging is 
not available), up to a total of 16, from CFHT U band in the near- 
ultra-violet to the 4 IRAC bands, but with the advantage of very 
deep and matched J, H and K data from UKIRT/UKIDSS DR5 
(c.f. Cirasuolo et al. in prep.). The appropriate corrections (equiva- 
lent to those used for COSMOS) are applied to obtain total fluxes 
in each band for the UDS field. Once all corrections are applied, 
the COSMOS and UDS magnitude/flux distributions are very well 
matched. 

Stellar masses are similarly derived for all galaxies within the 
underlying population at z = 0.84, assuming that they are all at 
z = 0.845 for simplicity (and that they are also representative of 
galaxies at z ~ 0.8); this thus provides a direct comparison with the 
masses estimated for the Ha population (after taking into account 
the completeness and contamination of the underlying sample - see 
Appendix A). 

Finally, the reader should note that even with a large number 
of photometric data points available, the stellar mass estimate of 
each individual source is found to be affected by a Icr error (from 
the multi-dimension x^ distribution) of ~ 0.15 dex, which results 
from degeneracies between the star formation time-scale r, age, 
extinction and, to a smaller extent, metallicity. However, these are 
the typical errors found in other literature estimates at 2 ~ 1, and 
tests show that the results of this paper are robust against the uncer- 
tainties in the stellar mass estimates. 

Other stellar mass esti mates in COSMOS have been presented 
by Mobasher et al. ( 2007| rj However, these have been obtained 
using a much less sophisticated methodjand are insufficiently ac- 
curate to undertake the analysis presented in this paper. Although 
the mass estimates of Cirasuolo et al. in UDS have been obtained 
using the same fitting method, there are some differences in the 
adopted photometric corrections and the SED templates used, so 
they are re-calculated here to ensure complete consistency between 
UDS and COSMOS. The reader is referred to Appendix B for a 
direct comparison between the stellar masses obtained in this pa- 
per and those from |Mobasher et al.H2007} and Cirasuolo et al. (in 
prep.). 



20 



^ The uncertainties in the luminosity distance (and the location of main 
spectral features) due to the filter profile are random and lead to much 
smaller uncertainties (^ 2 per cent) than other errors and degeneracies in 

the models. 

^ Improved mass estimates have been obtained by |llbert et al.|J2010| using 
the similar spectral energy distribution (SED) fitting method used here, but 
these are not publicly available. 

^ These were obtained from rest-frame V luminosities, My , and their ex- 
pected relation with stellar mass using colours drawn from the best spectral- 
type fit (c.f.|Mobasher et al.|2007). 



10 



o 

Q. 



2 



A*' 

X-'V Field 




Rich 
Groups/ 
Clusters 



Groups 



2 5 10 20 50 100 

Ec (Mpc-2) at 2 = 0.84 

Figure 1. The empirical relation between the correlation length, rp, and lo- 
cal density (corrected local projected density), derived for the Ha emitters. 
This provides an intuitive interpretation of both local densities and rg, giv- 
ing a better insight into the nature of the regions being probed and studied 
by the survey. This is found to be well fitted by a simple linear fit, given 
in Section 2.7. Bin sizes are chosen to guarantee that sub-samples have the 
same number of galaxies. 



2.6 Local density estimates 

Local densities for each source (Ha emitters and galaxies in the 
underlying population samples) are estimated by conducting a 10th 
nearest neighbour analysis. Briefly, for each source, the angular dis- 
tances to all other nearby sources are computed. These are then 
used to determine the projected distance to the lOth-nearest neigh- 
bour (in Mpc, using a scale of 7.64 kpc/arcsec at 2 = 0.845) and 
hence to estimate the local projected surface density. Sources near 
the edges of the fields and/or near masked regions will naturally 
be biased towards lower local density estimates (~5 per cent of 
them are in this situation). Where possible, this is fully tackled by 
using high-quality photo-z catalogues which exist well beyond the 
HiZELS NBj coverage. For regions where that is not possible (~ 4 
per cent of the UDS emitters are found near the edges of the photo-z 
data-set or large masked regions, and < 1 per cent of the COSMOS 
galaxies lie near masked regions), a simple correction is applied 
based on the fraction of area lost to the 1 0th nearest neighbour. 
Tests show that this does not introduce any bias. 

Local density estimates derived using the underlying (photo- 
z) sample are corrected by assuming the contamination and com- 
pleteness fractions derived in Appendix A. In general, given a con- 
tamination fraction, Co, a completeness fraction C, and a measured 
projected surface density E, the corrected local surface density is 
estimated as: 

Ee = Exi^. (1) 



C 



When Equation 1 is used to correct densities estimated using dif- 
ferent background populations (selected on the basis of photomet- 
ric redshifts - see Section|Z2l the different estimations are found to 
agree very well, albeit that the scatter increases for wider (in photo- 
z cut) underlying population samples (as contamination becomes 
increasingly dominant. However, the scatter between density mea- 
surements obtained using different photo-z samples, or derived us- 
ing different density estimators (e.g. 5th versus lOth nearest neigh- 
bour), is always significantly smaller ( <, 0.18 dex) than the bin 
sizes used in this work (~ 0.30 dex). Thus, as Appendix A demon- 
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Figure 2. The on-sky distribution of tlie star-forming Ha emitters found at z = 0.84 in the UDS (left panel) and COSMOS (right panel) fields when compared 
with the local projected density field which can be linked with qualitative environments - c.f. Section [Tt] Sources classed as likely and potential AGN in 
Gam et al. (2010) are not shown. Circles mark the position of extended X-ray emission from confirmed groups and clusters (from [Knoguenov et al.|2007] 
[2010 1 within the narrow-band redshift range, and are scaled to reflect the measured X-ray luminosity (in a log scale). Some lower X-ray luminosity groups are 
identified with black circles, while richer X-ray groups are plotted in white for good contrast. Note the very rich structure in the COSMOS field (see Section 
2.8), providing a unique opportunity to probe the densest environments. 



strates, the main results of the paper do not depend on the choice 
of the underlying population; the samples used for the analysis are 
the ones presented in Section[2^ 



Figure [T] shows how ro correlates with Ec. The trend is found 
to be very well approximated (forro inh^^ Mpc and Ec inMpc^^; 
see Figure[T|l by ro ~ 0.4 x Ec + 1.6. 



2.7 The clustering-density relation and the physical 
interpretation 

Estimating local surface densities can quantitatively distinguish 
higher density regions from lower density regions: but how do the 
derived surface densities relate to qualitative environments such as 
fields, groups or even clusters? Here, an attempt to associate qual- 
itative environments with ranges of Ec will be presented, by using 
the real space correlation length, ro. 

Following the procedures fully detailed in Sobral et al. 
(2010a), the entire sample of Ha emittersjis split based on cor- 
rected local surface density, Ec, in order to compute the correlation 
lengtnj of each sub-sample (using the exact relation between the 
angular and spatial correlation functions; Sobral et al. 2010a). 



° The sample of Hq emitters is used instead of the photometric-redshift- 
selected population as the real-space correlation length can be well esti- 
mated (as detailed in Sobral et al. 2010a) without major errors being intro- 
duced by the uncertainty in the exact redshift distribution and contamina- 
tion. 

'^ The analysis was carried out fixing 7 = —1.8, since individual sub- 
samples are too small to accurately determine it. This value of 7 provides 
a good fit to the overall sample (though note that, as discussed in Sobral 
et al. (2010), for narrow-band surveys this does not necessary translate to 
/3 = —0.8 as the Limber approximation would imply, and that the ex- 
act transformation between the angular and the 3-D correlation function is 
used). 



Having been able to find an empirical link between Ec and 
the correlation length of the Hq emitters, one can associate dif- 
ferent qualitative environments with local surface density (based 
on ro). Studies such as Yang et al. (20051 found that galaxies in 
the poorest group environments typically have ro ~ 6h^^Mpc, 
while the rich groups have ro « 15h~^Mpc; for cluster envi- 
ronments, studies such as Mo et al^ ( [1996 ) find typical correla- 
tion lengths of ro > 15 h~ Mpc. For the purposes of the analysis, 
three main environments will be considered: field (including poor 
and rich field), poor groups and rich groups/clusters. Field envi- 
ronments are defined to be those with Ec < lOMpc"^ as these 
correspond to ro < 6h~^ Mpc. Poor group environments are con- 
sidered to be those for which Ec = 10 — SOMpc^'^, as samples 
with these densities yield the range of ro found for such groups 
(6 < ro < 15 h^^ Mpc). Rich groups/clusters are considered to be 
those at even higher densities, corresponding to ro > 15 h~^ Mpc, 
for which Ec > 30Mpc"^. 



Of course, these simple classifications do not attempt to be 
fully accurate, nor do they aim to identify the exact transitions be- 
tween what should be considered field, groups or cluster environ- 
ments on their own. Nevertheless, they provide a valuable tool to 
associate the estimated local densities with different qualitative en- 
vironments. 
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Figure 3. Left: The dust-extinction corrected SFR (Ha, SFRs obtained after applying tlie dust extinction from Garn et al. 2010) as a function of stellar mass 
for both the original S09a Ha sample (EW > 50 A) and the possible Ho emitters with EW < 50 A, showing that the median SFR is uncorrected with stellar 
mass. Error bars present the 1 a Poissonian errors (ztcrm/vN). Right: the specific SFR (sSFR; SFR divided by stellar mass) falls with stellar mass, for both 
samples, although the steepness of such fall at the lowest masses is at least partially driven by the flux limit of the survey (dashed line; see discussion in main 
text). Both panels show very clearly how the nature of the samples with different EWs differs, with the S09a sample being made of galaxies with lower stellar 
masses, higher SFRs and thus much higher sSFRs at all masses, while the lower EW sample is mostly made of massive galaxies, presenting low sSFRs. 



2.8 Confirmed clusters and groups in COSMOS and UDS 

The on-sky distribution of Ha emitters can be seen in Figure[2] to- 
gether with the different densities in which they reside in both the 
UDS (left panel) and the COSMOS (right panel) fields. By prob- 
ing a very wide area, the Hq survey is able to target a relatively 
wide range of densities. Particularly, by chance, the survey is able 
to probe a luminous and massive (M500 ~ 10^*'' Mq) X-ray se- 
lected cluster in COSMOS at z = 0.835 (see Figure[2](, along with 
other nearby lower luminosity X-ray groups at 2: = 0.85 which 
have already been reported by the COSMOS team ( |Finoguenov| 
[et al.|2007[|Lea uthaud et al. 2010). These confirmed group/cluster 
regions are identified as high local projected density regions, with 
the calibration from §2.7 clearly classifying them as group/cluster 
regions: this not only confirms that this study is able to probe a 
very wide range of environments, from the poorest fields up to 
rich groups/clusters, but it also backs-up the qualitative environ- 
ment calibration used. 

[Finoguenov et al.] ( [2010| ) reports on a search for groups and 
clusters in the UDS field; their catalogue is used to confirm a rea- 
sonably rich X-ray group in the UDS field, along with 2 other lower 
luminosity selected groups within the narrow-band redshift cover- 
age of Ha (see Figure [2]l. This further confirms the wide range of 
densities probed by the survey, not only in the COSMOS field, but 
also in the UDS field. 



3 MASS-DEPENDENCES AND DOWNSIZING 

The complete sample of Ha emitters at z — 0.84 present typical 
stellar masses of M= 10^"'^ M©, while the general population of 
galaxies at the same redshift are redder in the rest-frame U-V, and 
present a wider range of stellar masses. 

As the left panel of FigurelSlshows, the median Ha star forma- 
tion rate (corrected for dust extinction following Garn et al. 2010) 
for the entire sample does not correlate with stellar mass. The lat- 
ter also reveals that the lower EW Ha emitters are typically much 
more massive and less active. This is further confirmed by com- 



puting the specific star formation rates (sSFR, SFR per unit stellar 
mass) of the entire sample: as the right panel of Figure [3] reveals, 
there is a significant decrease in sSFR with stellar mass for the two 
samples of Ha emitters, with the lower EW Ha emitters present- 
ing the lowest sSFRs (stressing the importance of including these, 
since they allow for a much wider range of masses and sSFRs to be 
studied). 

The observed sSFR-mass trend is similar to the one found by 
studies such as |Feulner et al.| pOOSj, but is contrary to other de- 
terminations which suggest a much flatter relation between mass 
and sSFR (e.g. "Dadd i et al.|2007l [Dunne et al.][2009l [Rodighiero[ 
|et al.|2010j ). One reason for this relates to a selection bias of the 
current sample due to the fixed SFR limit (dashed line in Figure 
[51: at the lowest masses probed, low sSFR galaxies fall below the 
Ha selection limit, biasing upwards the median sSFR determined. 
To estimate the mass below which this bias becomes significant, 
the survey limit of SFR > 3 M0 yr^^ can be combined either with 
the typical median sSFR~ 10~^ yr~^ argued in the above papers, 
or with the 6 Gyr age of the Universe at 2; ~ 0.8, by consider- 
ing galaxies which form stars at a fairly constant rate across their 
cosmic history. These both suggest that a detection bias will set in 
below ~ 10^" M0 (a value confirmed by the analysis later in this 
section). At lower masses, the galaxies that are detected are mostly 
those with unusually high sSFRs (i.e., bursty). 

This selection bias will therefore explain some or all of the 
steeper slope at low masses in Figure 3, but it cannot explain the 
strong trend at higher masses. Here, the difference from (and be- 
tween) previous literature determinations largely arises from the 
definition of "star-forming galaxies". In the current work, all galax- 
ies with Ha detections are classified as star-forming (SF) galaxies, 
leading to a clean selection cut in SFR. In contrast, previous works 
have often applied additional criteria to select the SF population. 
For example, Daddi et al. analyse 24 [im and ultra-violet data (i.e. 
broadly a SFR-selection, as in the current work) for a population 
of BzK-selected galaxies, but then they remove from considera- 
tion the galaxies classified as 'passive BzKs' (a simple colour cut). 
This colour cut broadly corresponds to a cut in sSFR of their galaxy 
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Figure 4. The fraction of galaxies forming stars with SFRs > 3 M© yr~ ^ 
when compared to the underlying population as a function of stellar mass 
for the entire sample, separately for the COSMOS and UDS fields, and for 
the S09a sample only. The figure shows that the fraction of star-forming 
galaxies with SFRs above the survey limit (SFR > 3 Mq yr^^) is Ri 40 per 
cent at ~ 10^" Mq, but falls sharply to as a function of stellar inasses. 
The EW cut of S09a results in a sharper decline of the star- forming fraction, 
mostly because it misses the most massive galaxies with low specific star 
formation rates. Errors are Poissonian, but they also include a 20 per cent 
contribution from the containination and completeness corrections applied. 



population, and has the effect of excluding the high-mass low-sSFR 
galaxies from their sample: it is not surprising that this then leads to 
a relatively flat mass-sSFR relation. Thus, the nature of the "mass- 
sSFR relation for star-forming galaxies", which has been widely 
debated in the literature, is critically dependent on exactly which 
galaxies one considers to be star-forming. 

By comparing the sample of H« emitters (star-forming popu- 
lation) with the underlying galaxy population it is possible to de- 
termine the fraction of star-forming galaxies as a function of stellar 
mass, which can provide an important insight into where most of 
the star formation is occurring dX z ^ 1. Figure H] presents the re- 
sults (after correcting the underlying population for contamination 
and completeness as a function of mass; see Appendix A). It reveals 
that while for M ~ 10^" M© the fraction of star-forming galaxies is 
~ 40 per cent, it continuously drops as a function of increasing stel- 
lar mass, effectively reaching per cent just below 10^^'^ Mq; the 
same results are independently recovered for COSMOS and UDS. 
Using the S09a sample produces a sharper decline as a function 
of mass, mostly because the EW selection misses massive galaxies 
with very low sSFRs. 

Figure PI also reveals an apparent drop in the fraction of star- 
forming galaxies at lower masses (< 10^'^Mq). This is simply 
a result of the SFR limit in the Ha sample that was discussed 
abovaj Nevertheless, specific star formation rate cuts can provide 
a view which broadly avoids this bias at the low mass end. Figure 
[5] presents the fraction of star-forming galaxies presenting sSFR 
above a certain value (points are only plotted if they are above the 
SFR limit of the survey, which avoids biases from the flux limit). 
The results reveal that the trend of a decreasing star-forming frac- 
tion with increasing stellar mass is recovered for galaxies with a 



* In order to properly probe star-forming galaxies with lower masses one 
requires a much deeper survey such as ROLES dGilbank et al. 12010). 



O logsSFR> -10.25 

■ log sSFR> -9.76 

^ log sSFR> -9.25 

j^ logsSFR > —8.75 

• log sSFR> -8.25 




10.0 10.5 11.0 
log Stellar Mass (M^) 

Figure 5. The fraction of star-forming galaxies as a function of stellar mass 
with a specific star formation rate greater than a given cut-off (corrected 
for completeness and contamination as a function of mass). This shows 
that the downsizing trend presented in Figure l4] is recovered regardless of 
the sSFR cut-off used, and allows for lower masses to be probed. En'ors 
are Poissonian and take into account the contamination and completeness 
from the photometric redshift selected population. Points are only plotted if 
above the SFR> 3 Mq yr~ ^ survey limit. 



wide range of specific star formation rates, implying that the re- 
sults are not simply biased by more passive-like, massive galaxies, 
dominating the very massive end of the star-forming fraction-mass 
relation. This analysis is also able to probe to lower masses, where 
the same trend is recovered. Nevertheless, it should be noted that 
the trends recovered for different sSFR cuts are not simple linear 
functions and, particularly for high sSFR cuts, the fraction of star- 
forming galaxies above a certain stellar mass is effectively zero up 
to the highest masses probed. 

These results indicate that mass-downsizing is already in place 
at z ~ 1, with a sharp decrease of the star-forming fraction as a 
function of stellar mass for the highest masses. Therefore, the most 
massive galaxies at 2 ~ 1 are mostly non-starforming, and even 
those that are forming stars present very low sSFRs. 



4 THE ENVIRONMENT OF Ha EMITTERS 

In this Section, the star-forming fraction is studied as a function of 
corrected surface density at z — 0.84. The main result can be seen 
in the left panel of Figure |6] which reveals that the star-forming 
fraction is relatively flat with density up to Ec ~ lOMpc^^ (with 
hints of a slight rise, but only at a ~ 1.5a level), but then 
falls sharply for much higher densities. This seems to be associ- 
ated with the qualitative transition between environments: the star- 
forming fraction (for galaxies producing stars at a rate higher than 
3 Mq yr^^) decreases with increasing density once group regimes 
are reached, approaching zero for the highest densities; this is in 
good agreement with the results presented by |Koyama et al.| l |2010[ > 
at a similar redshift. The trend is the same in both COSMOS and 
UDS fields (even though the UDS field is not able to probe the high- 
est densities found in the COSMOS X-ray cluster), revealing that 
this is not caused by a bias in one of the fields, nor due to the se- 
lection of the underlying population (this is very unlikely since the 
COSMOS and UDS photometric redshift catalogues are completely 
independent, obtained with different codes, and by different teams). 
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Figure 6. The fraction of the H«-emitting galaxies when compared to the underlying population at the same redshift (left panel) and the median Ha star 
formation rate of Ha emitters (right panel), as a function of local projected galaxy surface density (for Ha emitters with SFRs > 3Mq yr~^). The error 
bars present the Poissonian uncertainty for each bin and in the left panel also include a 20 per cent contribution from the completeness and contamination 
connections applied. At the lowest densities, the star-forming fraction is relatively constant and only increases slightly, from ~ 30 per cent at the lowest field 
densities probed to ~ 35 per cent across the field environment. However, for higher densities, Sc > lOMpc^^, there is a steep decline of the star-forming 
fraction down to the highest (rich groups/cluster) densities probed, where the star-forming fraction is only Ri 5 per cent. The results also show that the typical 
(median) star formation rate of Ha emitters increases continuously from the lowest densities to Sc ~ 50 Mpc~^, although it appears to fall again in the very 
richest environments. This corresponds to emitters residing in groups and cluster outskirts presenting typical SFRs which are about 2 times higher than those 
residing in the field. 



It should be noted that the narrow-band filter is not a perfect top-hat 
and that some group/cluster members may have their Ha line be- 
ing detected towards the wings of the filter, affecting the likelihood 
of detecting them as star-forming. However, particularly because 
two independent fields are considered in the analysis, this effect is 
likely to even out in the statistics, with different groups detected at 
all points on the filter profile. Furthermore, as detailed in Appendix 
A, the trend is recovered regardless of the choice of underlying 
populations: only small differences are found once one applies the 
completeness & contamination corrections given by Equation 1. 

The right panel of Figure[6]shows that the median star forma- 
tion rate of Ha emitters rises with increasing local density, at least 
up to group densities, although there are hints that this relation does 
not hold for highest (rich groups/cluster) densities, consistent with 
previous results which find that, for rich environments, the typi- 
cal SFR decreases with S (e.g. |Poggianti et al.|2008i l. The results 
in this paper are also consistent with Sobral et al. (2010), which 
found that the correlation length of star-forming galaxies rises with 
Hq luminosity. The relation between the median SFR and density 
is seen in both COSMOS and UDS, and is also observed for differ- 
ent Ha luminosity limits adopted or when the analysis is restricted 
to spectroscopically confirmed Ha emitters. 

Moreover, as Figure [7] shows, the connections between envi- 
ronment and star formation activity are likely to be related with the 
star formation mode: while star formation in the field is dominated 
by non-interacting galaxies, at rich group and cluster densities, the 
star formation seen at z = 0.84 in Ha is dominated by potential 
mergers; with this correlation being found at a 2Aa level. In fact. 
Ha star-forming galaxies found in the X-ray core of the massive 
cluster in COSMOS are all mergers. 

The environment seems to be connected with a change in the 
star-forming fraction, the typical star formation rate and the nature 
of Ha emitters at 2: = 0.84; how are these changes manifested 
in the Ha luminosity function? While it is nearly impossible to 
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Figure 7. The potential merger and non-merger fractions as a function of 
projected local surface density (errors are Poissonian). This reveals that the 
fraction of potential mergers is reasonably flat at ~ 33 per cent for field 
densities but is then found to increase continuously with density, reaching 
~ 100 per cent at rich groups/cluster densities. 



accurately determine the volume occupied by galaxies residing in 
different environments, one can arbitrarily normalise any obtained 
Ha luminosity function to investigate changes in L'^^ and the faint- 
end slope, a: this approach can provide the insight that one needs 
to further investigate the relation between star formation and envi- 
ronment. 

Figure Is] presents the normalised Ha luminosity functions^ 

^ Ha luminosity functions are computed using the entire sample of Ha 
emitters, which includes the low EW sample. Conclusions remain un- 
changed even if the low EW sample is excluded. 
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Figure 8. The normalised Ha luminosity function for different environ- 
ments, with the best-fit Schechter functions. While the characteristic lu- 
minosities, L^^, of the different luminosity functions are in very good 
agreement within the errors, the results show a significant difference in the 
faint-end slope, a, as a function of local surface density. As shown in the 
inset, this is found to be very steep {a Si —1.9) for the lowest densities 
(poor field, E < 5Mpc~'^), shallower for medium densities (rich field, 
5 < E < lOMpc"^, with a fa —1.5) and very shallow for gr oups and 
clusters (S > lOMpc"^), with a Ri —1.1. Recent results from |Koyama| 
|et al.H 2010), which studied one cluster and its outskirts at 2 = 0.8 (there- 
fore equivalent to group and cluster densities) fully agree with the Ha lu- 
minosity function derived for the Risame local projected densities. 



appropriately weigtited as discussed in Section 2.3) for different 
local projected densities. By fixing </>*, both LJ^^ and the faint- 
end slope, a, are allowed to vary in order to obtain the best 
Schechter function fit to each luminosity function. The results re- 
veal that while the value of LJj^ is consistent within the errors 
for all environments probed by the HiZELS survey at z = 0.84 
(log Lhc = 42.47 ± 0.05, 42.44 ± 0.07, 42.45 ± 0.09 erg s'^ for 
poor fields, rich fields and groups and clusters, respectively), there 
is a strong difference in the steepness of the faint-end (see the inset 
of Figure [HI. The faint-end of the LF becomes continuously shal- 
lower with increasing density, from a — — 1.9±0.2, for the lowest 
densities, to a = —1.5±0. 2 for rich fields and a ~ — l.l±0.2for 
groups and clusters. This clearly implies that the enhancement in 
the typical SFRs with environment directly results from the strong 
depletion of the faint-end and the shallow a in rich environments 
relative to the field, rather than an increase in the bright-end. This 
also reveals that the steepness of the faint-end of the Ha LF in 
S09a is completely driven by field galaxies. The reader is referred 
to Section 6.2 for a discussion and possible interpretations of these 
results. 

Figure [8] also provides a direct comparison with results from 
|Koyama et al.| ( (2010| >; these authors studied a high density region 
at z ~ 0.8, probing densities which are directly comparable with 
groups and cluster densities. The data are corrected for extinction 
using the relation presented in Garn et al. (2010) for a direct com- 
parison. Koyama et al. found their data to be reasonably consistent 
within the errors with the total S09a LF (which contained all envi- 
ronments), but their data fully support the difference in the faint- 
end of the luminosity function between field and group/cluster en- 
vironments found here, and completely agrees with the results de- 
rived for the same densities, supporting the results which show that 
the environment seems to be responsible for setting the faint-end 
slope of the LF, at least at z — 0.84. A qualitatively similar result 



has been found for the 0024 super-cluster by Kodama et al. (2004) 
at 2 = 0.4; however, while studies at even lower redhift (see e.g. 
|Balogh et al.|2002| > have found different and mostly discrepant re- 
sults regarding a. It is clear that a comprehensive and large local 
study is still needed to really detail the role of the environment on 
the faint-end of the Ha LF in the local Universe, to provide the 
low-redshift benchmark that can be used to look for any cosmic 
evolution.. 



5 THE NATURE OF STAR FORMATION AND THE 
MASS-ENVIRONMENT VIEW 

5.1 The star-forming fraction in the 2-D mass-environment 
space 

It has been shown that star formation has strong dependences on 
both mass and environment at z ^ 1. Nevertheless, with environ- 
ment and stellar mass expected to be somewhat correlated, how 
related are the dependences of star formation on mass and environ- 
ment and are these intrinsically connected? Figure [9] presents the 
mass-environment 2-D space for Ha emitters and for the entire un- 
derlying population. This reveals the expected trend of an increase 
in the typical stellar mass with increasing local projected density 
for both populations of galaxies (with a change of ~ 0.5 dex in the 
median mass from low to high densities). 

To attempt to distinguish any potential inter-dependence of the 
roles of mass and environment on star formation at z ~ 1, the 
analysis presented in the previous 2 sections is repeated in the 2-D 
space of mass and Ec. The right panel of Figure[9]presents the star- 
forming fraction evaluated within this space, in regions where this 
can be calculated (see caption for details). While this clearly illus- 
trates the overall behaviour which has already been presented, the 
fine details suggest significant connections between the dependence 
of the star- forming fraction on mass and environment. In particular, 
the results show that the weak increase of the star-forming fraction 
with density within the field regime and its subsequent sharp de- 
cline at higher densities is predominantly driven by galaxies of typ- 
ical mass (~ 10^" — 10^"'^ Mq) in the sample (please note that, as 
described in Section 3, there is a significant incompleteness effect 
at the lowest masses; this is indicated by the horizontal filled line 
in the Figures). Interestingly, whilst the mass-downsizing trend is 
driven by galaxies in the dominant field and group environments, 
the results reveal an increase in the fraction of star-forming galax- 
ies with density at the highest masses and, equivalently, an apparent 
reversal of the mass-downsizing trend in the richest enviroimients. 

Earlier, it was shown that the fraction of potential mergers in- 
creases with local density. To investigate whether this change in 
the potential nature of star formation can shed some light onto un- 
derstanding what is happening at relatively high masses and high 
densities, the left panel of Figure [To] presents how the fraction of 
star-forming galaxies that are classified as potential mergers (in the 
COSMOS field) varies within the 2D mass-environment space. It 
clearly reveals a strong increase in the fraction of potential merg- 
ers with environment, and, particularly, that star formation at the 
highest densities is mostly driven by potential mergers, regardless 
of mass. It is therefore clear that the behaviour of the star-forming 
fraction within the studied 2D space reflects the role of potentially 
merger-driven star formation for very high densities, where merg- 
ing events are more likely to happen. This is shown more clearly 
in the right panel of Figure [TO] which repeats the analysis of Fig- 
ure 9, but now restricted to isolated star-forming galaxies only (i.e. 
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Figure 9. Left: The distribution of both Ha emitters (filled squares: emitters with EW>50A; filled triangles: emitters with EW<50A) and galaxies 
within the underlying population (dots) in the 2D mass-environment plane. Vertical lines represent the qualitative transitions between field, groups and rich 
groups/cluster environments, while the filled horizontal line indicates the approximate mass for which the star-forming population becomes significantly 
incomplete due to its star formation rate limit (see Section 3). This reveals the expected correlation between median stellar mass and local density (larger 
symbols). Note that the median of the underlying population is evaluated over all galaxies, while the star-forming population is incomplete below the horizontal 
line, hence the apparently higher median mass for star-forming galaxies. Right: The fraction of galaxies that are forming stars as a function of both stellar 
mass and environment. This is calculated in regions of the 2-D space for which it is possible to gather 5 or more background galaxies in a (Gaussian) smoothing 
radius up to Ic of the typical errors in mass and environment - regions with less than 5 background galaxies are not shown. The latter reveals that the overall 
behaviour of the star-forming fraction as a function of stellar mass or environment depends, to some extent, on the environment or stellar mass, respectively. 
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Figure 10. Left: The fraction of potential mergers within the star-forming population of galaxies (in COSMOS). This reveals that star formation occurring 
at the densest environments is potentially merger-driven. Right: The non-merger driven star-forming fraction, revealing that the trends with environment and 
mass become even stronger once potential mergers are removed from the sample. This is calculated in regions of the 2-D space (COSMOS only) for which it 
is possible to gather 5 or more background galaxies in a (Gaussian) smoothing radius up to la of the typical errors in mass and environment - regions with 
less than 5 background galaxies are not shown. 



excluding all potential mergers). It is found that the fraction of non- 
merging star-forming galaxies peaks roughly at the mass complete- 
ness limit (for the SFR limit of the current sample, ~ 10^* Mq), 
and in the transition between fields and groups, but shows strong 
general declining trends with both increasing mass for all environ- 
ments and with increasing local density for all masses. This reveals 
that non-merger driven star formation is strongly suppressed in both 
rich groups/cluster environments and for high stellar masses, im- 



plying that once potential mergers are neglected, stellar mass and 
environment both play separate and important roles. 



5.2 The mass-SFR-environment dependence 

Section 4 showed that the typical SFR increased with density 
at least up to groups/cluster outskirts densities. Here this is in- 
vestigated further, by examining how the median star-formation 
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Figure 11. The dependence of the median SFR on environment (Sc) for 
star-forming galaxies with different stellar masses. Error bars present the 
1 a Poissonian uncertainty (it(Tm/'\/N). The results show that while the 
increase of the median Hq SFR with environment is valid for galaxies with 
stellar masses in the range Si ]^q9. 0-10.6 Mq^ the most massive galaxies 
in the sample (stellar masses > 10^*^'® M©) do not have their typical SFR 
enhanced with environment. 



rates vary with environment for sub-samples of star-forming galax- 
ies with different stellar masses. Figure 11 reveals that while a 
strong SFR-environment trend is seen in low to moderate stel- 
lar mass galaxies, it does not hold for the most massive galaxies 
(M > 10^"'® Mq), for which the median SFR seems to remain un- 
changed with increasing local density. Since the typical stellar mass 
of galaxies increases with density, in terms of sSFRs one finds that 
for the most massive galaxies, the median sSFR decreases with in- 
creasing local density (at a 2a significance level), whilst this trend 
is continually reversed with decreasing stellar mass, becoming a 
significant positive correlation (3.5cr) at the lowest masses probed. 
It should be noted that although the most pronounced difference 
in the SFR-E behaviour is in the densest environments, where the 
merger fraction is high, the same trends are recovered (albeit weak- 
ened by ~ 0.6it) if the analysis is conducted using only non-merger 
star-forming galaxies. Therefore, merging activity by itself is not 
enough to explain the enhancement of SFR in low-to-medium mass 
star-forming galaxies. These results imply that there is some inter- 
dependence in the influence of mass and environment on galaxy 
evolution. 



5.3 The colours of star-forming galaxies in the 2-D 
mass-environment space 

The nature of the Ha emitters occupying different regions in the 
2D space being probed can be made even clearer by looking at the 
colours of emitters. For z ~ 0.8, this can be simply done using 
{R — z) because these span the 4000 A break, and are sufficiently 
close so that dust extinction is unlikely to be a significant problem. 
Figure [12] (left panel) presents how the median colours of the Ha 
emitters vary across environment for different stellar masses. The 
results suggest that stellar mass is the most important quantity for 
determining the 4000 A break colour (and whether the Hq emitter 
is consistent with already being in the red sequence or still occupy- 
ing the blue cloud) of star- forming galaxies. Figure [T2] also reveals 
that the environment plays a much weaker role for determining the 
colour of star-forming galaxies, even though the reddest Hq emit- 



ters are found in group environments (but those are also very mas- 
sive). The 2-D colour distribution is found to be qualitatively the 
same for both potential mergers and non-merging galaxies, show- 
ing the same trends, although potential mergers are consistently 
bluer (by ~ 0.2 in colour) than the non-mergers throughout the 
2-D space studied. 

Further information can be obtained by comparing the colour 
distribution of the star-forming galaxies with that of the underlying 
population. As the right panel of Figure[T2]shows, the star-forming 
population is bluer than the underlying population at all masses 
and environments. However, the results indicate that even though 
the main colour predictor for star-forming galaxies is stellar mass, 
the difference in colour between the star-forming population and 
the underlying population of galaxies depends on both mass and 
environment. Particularly, the colour difference is greatest at lower 
masses (above the mass completion limit - see Section 3), where 
the underlying population contains a mix of blue and red galax- 
ies, while the star-forming population has very few red ones, thus 
leading to the difference. However, such difference is negligible at 
the very highest masses, where even star-forming galaxies are very 
red. Moreover, the fraction of redder galaxies at moderate masses is 
found to be higher in denser environments for the underlying popu- 
lation (the underlying galaxy population has a weak trend of colour 
with increasing density), but the star-forming population has essen- 
tially no trend with environment for a fixed mass (see left panel of 
Figure [T2|, thus resulting in the environmental trend of the colour 
difference. 



6 DISCUSSION 

It is clear that the bulk of the most massive galaxies that are already 
in place at z = 0.84 are non-star-forming (star formation rates 
lower than the HiZELS limit, 3 M0 yr~ ^ ) and many of those which 
are star- forming have low sSFRs and have stronger 4000 A breaks 
(are redder) than lower mass galaxies. These results clearly point 
towards mass downsizing already being in place for the general 
population of star- forming galaxies studied at z ~ 1, although po- 
tential merging activity (which becomes more likely at higher den- 
sities) can be responsible for triggering some star formation even 
in massive, apparently red-sequence galaxies. 

In addition to this, the environment has a significant role in 
the nature and evolution of star-forming galaxies at 2 ~ 1. For lo- 
cal surface densities associated with the field regime, and for mod- 
erate stellar masses, the fraction of star-forming galaxies is rela- 
tively flat with increasing envirorunental density, while the median 
star formation rate of Ha emitters is found to increase with den- 
sity. However, for higher densities, the star-forming fraction falls 
sharply towards ~ per cent for the highest rich group/cluster den- 
sities. The environment is also found to play an important role in 
setting the faint-end slope of the Ha luminosity function, with the 
results showing that under-dense regions have Ha LFs with very 
steep faint-end slopes, while the highest densities have a shallow 
faint-end slope - this confirms analytic predictions such as those 
presented by |Trenti & StiaveUi| j2008^ . If this holds at different 
epochs (particularly, at higher redshift), then it suggests that the 
steepening of the faint-end slope with redshift that has been argued 
by many authors (for both Ha and UV selected samples, e.g. [Hop-] 
|kins et al.|2000l|Bouwens et al.|2007|[Reddy & Steidel|2009^ , may 
be mostly a consequence of structure formation, since at high red- 
shift groups and clusters of galaxies become rarer and rarer and the 
field environment dominates. At lower redshift, on the other hand. 
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Figure 12. Left: The distribution of tlie median R-z colour (probing the 4000A break at 2: = 0.84) within the mass-density 2D space for the Ho star- 
forming galaxies, revealing that stellar mass is the main colour predictor as the environment only correlates weakly with colour. Right: The difference in the 
median R-z colour between the population of star-forming galaxies and the complete underlying population, revealing that although star-forming galaxies are 
consistently bluer than the underlying population, the colour difference is a clear function of mass and also depends on the environment. 



the discrepancies in the faint-end slope discussed in S09a are likely 
to be caused by cosmic variance, with different surveys targeting 
regions with different densities. 

When one combines the dependences of star formation on 
mass and environment, a sharper picture arises: this shows that the 
increase in the star-forming fraction in the field regime is dom- 
inated by median mass (~ 10^" — 10^'^'^ Mq), essentially non- 
merging galaxies, and that the bulk of star-forming galaxies in 
high density environments are red, massive galaxies undergoing a 
merger episode. It is also found that the fraction of non-merging 
star-forming galaxies strongly declines in rich groups/cluster den- 
sities for all masses, and for high masses for all environments, 
implying that both mass and environment independently suppress 
non-merging star-formation. However, this also shows that the en- 
hancement of the typical SFRs with increasing density is only hap- 
pening at z ~ 1 for galaxies with masses lower than 10^"'' M0 
- for the most massive star-forming galaxies the environment does 
not seem to alter the typical SFRs and their properties seem to be 
mostly set by their stellar mass. 



6.1 Comparison with other studies 

Studies at z ~ 1 and above have found apparent discrepant re- 
sults for the relation between star formation and environment, with 
some authors claiming to have found an inversion in that relation 
when compared to the local Universe. In particular, |Elbaz et al.| 
( |2007 ) study the star formation rate-density relation at 2 ~ 1 and 
find that, instead of declining, SFRs of star-forming galaxies are 
enhanced with local projected density; [Cooper et al.| pOOS} find 
similar results. These studies typically probe field to poor group 
environments. On the other hand, |Poggianti et 31.^2008) presents a 
detailed study of a large sample of relatively high redshift groups 
and clusters (up to 2 ~ 1). Contrarily to [Elbaz eFal. ( 2007 1 and 
[Cooper et al"] ( (2bo8l l, |Poggianti et al.| ( (2bo8l l find that SFRs decrease 
with environment in the high density regime (but the results agree 
with a rise in the mean SFRs - which are [Oil] based - with density 
up to S ~ 30 Mpc"^, and with the fall after that). 



|Patel et ar]j2009[ > studied a 2 ~ 0.8 cluster and the surround- 
ing regions (using a mass limited sample, but selecting only rela- 
tively massive galaxies, and making use of MIPS 24 |j.m data) to 
show that the specific star formation rate decreases with increasing 
density at 2 ~ 1 and argued that star formation is suppressed with 
density even at that epoch. The authors comment on the big dis- 
crepancy of their results when compared to |Cooper et al.| ( (2008) or 
[Elbaz et ar| ( [2007[ >, and, whilst they argue that the differences could 
be arising from a combination of different effects, they are not able 
to reach a definitive answer. [Feruglio et al.[ pOTOb also study the 
role of the environment on star- formation and, using a sample of lu- 
minous and ultra-luminous infrared galaxies (LIRGs and ULIRGs) 
in COSMOS, find that for such massive galaxies there is not an 
increase in SFR with density up to 2 --^ 1. 

Fortunately, the large sample and the wide range of environ- 
ments and stellar masses within the HiZELS sample provides an 
ideal opportunity to attempt to reconcile these results. Indeed, one 
finds, just like [Elbaz et al.[f2007^ , that the environment enhances 
SFRs at least up to groups/cluster outskirts environments, but then 
decreases at higher densities as found by Poggianti et al. Further- 
more, if one looks at how the sSFR of massive galaxies (compa- 
rable to those studied by Patel et al. 2009) varies with environ- 
ment, the same trend is recovered: sSFRs of massive galaxies de- 
crease with E, mostly because median SFRs of massive galaxies 
are largely insensitive to environment, but the median mass in- 
creases at higher densities, leading to a lower median sSFR. Note 
that such a trend is reversed for lower mass galaxies, where the in- 
crease in SFR with density is fully valid; this reveals that the results 
from Patel et al. (a mass-selected sample), or Feruglio et al. (2010b 
(a sample of massive LIRGs and ULIRGs) are highly dependent 
on the stellar mass selection, and that these do not agree with re- 
sults from [Elbaz et al.[j200 7) mostly because galaxies in the later 
study are selected from their star formation emission lines (there- 
fore bluer, less massive galaxies). 

Furthermore, [Poggianti et aE] ( [2008^ also finds that the frac- 
tion of star-forming galaxies falls continuously for densities signifi- 
cantly larger than E ~ 10 Mpc~^ (the study is not sensitive to field 
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densities). Similar results are recovered by Koyama et al. (2010), 
who studied a cluster al z = 0.8 and found that the star-forming 
fraction (Ha emitters) falls from groups to the cluster core, while 
[Feruglio et"ar]|2010 ) finds that the fraction of LIRGs and ULIRGs 
also decreases with increasing density up to 2; ~ 1. These are the 
same trends found in this paper, but the current work extends this 
to much lower densities and, overall, probes a wider range of den- 
sities. 

Together, the results in this paper are able to simultaneously 
agree with previously apparently contradictory results, revealing 
that such discrepancies are caused by the studies focusing on galax- 
ies with significantly different masses and, most importantly, resid- 
ing in different environments. 

The analysis presented in this paper contributes towards im- 
proving our understanding of galaxy evolution and the key epochs 
for significant transitions. While the median SFR is enhanced with 
density up to moderate densities, for a star-forming selected pop- 
ulation, it has also been shown that the behaviour of that relation 
depends on stellar mass. Furthermore, while merging does become 
important at the highest densities, it is also found that the SFR-E 
relation is still recovered even when considering non-mergers only. 
This surely implies that the physical processes through which the 
environment enhances the typical SFRs cannot be fully explained 
by (major) mergers. Therefore, galaxy harassment, or the accelera- 
tion of the intergalactic gas infall in rich groups and cluster outskirts 
are likely to play a significant role in enhancing star-formation 
at high densities. Since these processes are more likely to affect 
galaxies with higher gas contents, it may be understandable that 
the trends are seen most strongly at low and medium masses. 

Furthermore, the star-forming fraction drops from fields to 
groups and clusters, revealing no qualitative evolution from the 
local Universe: very rich environments already have a low star- 
forming fraction even at z = 0.84 (and the star formation activity 
occurring at high densities is potentially merger-dominated). Par- 
ticularly when one excludes potential mergers, the environmental 
trend of a declining star-forming fraction for high densities is re- 
covered at all masses, revealing that the environment plays an im- 
portant role which is roughly independent from that of stellar mass. 
This is in good agreement with the results presented by |Cooper| 
|et al.| pOTO}, at a similar redshift, and qualitatively similar to the 
clear independence of the quenching roles of stellar mass and en- 
vironment found by peng et al.| ( |2010p using SDSS. Since HiZELS 
is able to agree with previous results and provide a much broader 
view, it further confirms that the densest regions at 2; ~ 1 are al- 
ready the sites where one finds predominantly red sequence galax- 
ies and low star-forming fractions, as found in previous cluster 
studies. However, with studies such as |Hayashi et al.|f2010 l show- 
ing that the star-forming fraction does not drop with increasing den- 
sity even for the densest regions at z ~ 1.5 (but that some of [Oil] 
emitters within the cluster core are already in the red sequence), it 
seems likely that the latter stages of the environment transforma- 
tion happen between z ~ 0.8, where one finds red Ha emitters at 
cluster cores, but a fraction of star-forming galaxies which already 
declines steeply, and 2 ~ 1.5, where such fraction seems to hold 
up to the highest densities; the same view seems to be recovered by 
using 24 urn data ( [Hilton etal.|2010^ . 



it is now clear that the bulk activity of such galaxies is occurring in 
the field and poor groups. 

Nevertheless, while the fraction of non-merger driven star- 
forming galaxies falls sharply at high densities, there is still some 
star formation found in the richest environments. The latter is 
merger dominated, and this is linked with the increase in the me- 
dian SFRs of Hq emitters and, to some extent, with the fall of the 
star-forming fraction with environment. 

Indeed, in group and cluster environments, ~ 80 per cent 
of star-forming galaxies are potential mergers. This is somewhat 
expected, as mergers are more likely to happen in dense envi- 
ronments. However, merger-driven star formation at z — 0.84 is 
likely to result in starbursts which quickly exhaust the available 
gas (which, for the most massive galaxies that are already reason- 
ably red, is likely to be low), and are expected to produce ellipti- 
cal galaxies. This would then further drive the decline of the over- 
all star-forming fraction with environment and would be consistent 
with a significant population of mergers at high densities at slightly 
higher redshifts, as argued by S09a to explain the evolution of the 
Ha luminosity function. Thus, as the densest regions increase the 
merger chances, they can lead to a much quicker migration of blue 
cloud galaxies into the red sequence. This would then establish a 
higher fraction of passive galaxies in the densest regions much ear- 
lier than in the field, where potential merging events are not domi- 
nant. 

At higher redshift (2 = 1.2), [Ideue et al.| ( [2009l l presented 
an environment analysis using [Oil] 3727 emitters selected from 
a narrow-band survey, accessing the potential role of interacting 
galaxies as a function of environment in the COSMOS field (but 
instead of morphologically classifying the galaxies, they looked at 
the fraction of isolated against non-isolated galaxies). The authors 
find an increase in the fraction of non-isolated galaxies with den- 
sity. However, while this is fully consistent with the results pre- 
sented in this paper and with the interpretation, it should be noted 
that |Ideuee t al. ( 2009 ) define non-isolated galaxies as galaxies sep- 
arated by distances of less than 80 kpc. Whilst this is not a robust 
criteria for identifying a merger, it also strongly correlates with the 
local density estimate, even if the merger fraction does not go up 
with density; therefore, an increase with density may be expected 
as a bias, regardless of the nature of the galaxies (real mergers or 
non-mergers). 

The results of this paper can also be used to interpret the influ- 
ence of the environment on the Ha luminosity function. In S09a 
the Ha luminosity function derived for just potentially merging 
sources was found to have a shallow faint-end slope, a. With this 
study it seems likely that the low a is a consequence of a higher 
density environment: galaxies in high densities will increase their 
merging chances and reside in potentially gas-rich regions at 2 ~ 1, 
making it hard to maintain low, isolated star formation rates. 

The results imply that higher densities enhance star formation 
in star-forming galaxies and that seems to be valid at all cosmic 
times; that's why denser regions have the most massive galaxies; 
however, because star formation has been enhanced, galaxies form 
stars quicker and sooner, and this leads to the mass downsizing 
observed since 2 = 1 down to the local Universe. 



6.2 The role of the star formation mode: potential mergers 
vs. non-mergers 

S09a find that non-merging galaxies (mostly disc galaxies) domi- 
nate the star formation rate density of the Universe at 2 ~ 0.8, and 



7 CONCLUSIONS 

The conclusions of this study can be summarised as follows: 

• The well-defined HiZELS sample of Ha emitters at z — 
0.845 ± 0.015 down to SFRs > 3M0yr~^ covering a large sky 
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area has been used to conduct a detailed study of the dependence 
of star formation on mass and environment. This is able to reconcile 
previous results and provide an integrated, clearer view of galaxy 
evolution at z ~ 1, particularly by taking advantage of the excel- 
lent underlying data in the COSMOS and UDS fields, but also be- 
cause the survey probes significantly massive, rich groups/clusters, 
providing a truly panoramic Ha survey over a very wide range of 
environments. 

• The star-forming fraction is found to strongly decline with 
stellar mass, from « 40 per cent for 10^° M© galaxies to effec- 
tively zero above 10^^'^ M©. Specific star- formation rates also de- 
cline continuously with stellar mass. Mass downsizing is therefore 
fully in place at z = 0.84. 

• The median SFRs of Ha emitters at z = 0.84 increases as a 
function of local surface density for both field and group environ- 
ments, but the trend is stopped for the highest densities, where it 
is likely to fall. This trend is driven by star-forming galaxies with 
stellar masses up to ~ 10^"'® Mq; higher-mass galaxies seem to 
have median SFRs which are mostly independent of the environ- 
ment they reside in. With this result being recovered even when one 
excludes potential merger-driven star-formation, it is likely that at 
least part of the environmental enhancement must come from other 
physical processes, such as galaxy harassment, or the acceleration 
of intergalactic gas infall in groups and cluster outskirts. 

• The fraction of galaxies that are forming stars is relatively flat 
with density within the field regime, but there is a clear, continuous 
decline once group densities are reached, resulting in the fall of the 
star-forming fraction from fields to rich groups/clusters, as seen in 
the nearby Universe and in studies focusing on the highest densities 
at 2 ~ 1. 

• The environment is responsible for shaping the faint-end of the 
Ha luminosity function: the faint-end slope, a is found to be very 
steep(a = — 1.9±0.2)forthepoorestregions,becomingshallower 
(a = 1.5 ±0.2) for rich fields and even shallower (a — — 1.1±0.2) 
for groups and clusters. 

• Star formation at the highest densities is dominated by poten- 
tial mergers, with likely implications on the nature and duration of 
the star formation activity, and suggesting that at least part of the 
environment quenching of star formation seen in the local Universe 
is driven by merging events at higher redshift. Neglecting potential 
merger-driven star formation results in stronger and mostly inde- 
pendent trends of a declining fraction of star-forming galaxies with 
both increasing mass and environment. 

• The 4000 A colour of star-forming galaxies correlates 
strongly with stellar mass, but only weakly with environmen- 
tal density. Massive star-forming galaxies residing in groups are 
mostly in the red sequence, while lower mass galaxies residing in 
the field occupy the blue cloud. 

These results show that both mass and environment play signifi- 
cant, inter-dependent roles in galaxy evolution, although such roles 
may become more independent if one neglects potential merger- 
driven star-forming galaxies. High density regions are linked with 
triggering high SFRs (for moderate mass galaxies), but this is likely 
to result in a quicker quenching of star formation as the gas is 
rapidly consumed, thus quickly transforming groups and cluster en- 
vironments into the predominant sites for gas-poor, passive galax- 
ies. On the other hand, the quenching processes linked with stellar 
mass (e.g. AGN) must also play an important role which is distinct 
from that of environment. It seems that the combination of high 
masses and high densities is the recipe for very effective and early 
galaxy formation, as even at z~ 1 the star-forming galaxies found 



with such combination are already on the red sequence, with very 
low sSFRs (even though most of them are potential mergers). 
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APPENDIX A: THE UNDERLYING POPULATIONS 

Al Completeness and contamination 

State-of-the-art photometric redshifts (Ilbert et al.'2009' Cirasuolo 
et al. in prep.) can be used to obtain samples of galaxies found at 
roughly the same redshift as the Ha emitters. In this Appendix, 
the completeness and contamination of the different photometric- 
redshift-selected samples are investigated and, after correcting for 
these, the robustness of the results presented in the paper is demon- 
strated. 

For COSMOS, where the photometric redshifts can be ob- 
tained more precisely (due to the availability of medium bands 
which can, for example, probe the 4000 A break extremely well), 
samples are defined selecting galaxies with 0.845 ± Az, with Az 
ranging from 0.005 to 0.065 (0.78 < Zphoto < 0.91) in steps of 
0.005 . Similar samples are defined for the UDS, but in steps of 0.01 . 
The samples are selected over the entire available area with the best 
photometric and spectroscopic data in both fields. 

For the COSMOS field, the zCOSMOS DR2 secure (> 99 per 
cent confidence) spectroscopic redshifts are used to study the com- 
pleteness (the number of galaxies selected by the photo-z cut that 
have spectroscopic redshifts within the top hat narrow-band red- 
shift range for the Ha emission line [0.83 < Zspec < 0.86] com- 
pared to the total number of galaxies with spectroscopic redshifts in 
that range over the same area) and the contamination (the fraction 
of galaxies in the photo-2; sample with spectroscopic redshifts that 
are confirmed to lie outside the narrow-band redshift range). The 
photo-z samples are all ~ 10 per cent spectroscopically complete 
and there are a total of 26(Lj spectroscopically confirmed galaxies 
in the full NBj filter range (over the entire COSMOS field). 

^^ From these, 182 are found over the NBj imaging coverage, with 
roughly half of these being detected by HiZELS as Ha emitters. 
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Table Al. A completeness and contamination study of the photometric red- 
shift selected samples obtained in the full COSMOS field using 2-COSMOS 
DR2 jLillyetafi2009 1. Completeness is estimated as the fraction of galax- 
ies wich have spectroscopic redshifts in the naiTow-band (top hat) filter 
range (0.83 < z < 0.86) which are recovered by the photo-2: cut (and 
after accounting for introducing all confirmed galaxies into each sample). 
The contamination is calculated by computing the fraction of galaxies with 
spectroscopic redshifts for which the redshifts are not within the nan'ow- 
band filter redshift range for Ho emitters at 2 = 0.84. Photometric-selected 
samples are all Ri 10 per cent spectroscopically complete. The samples 
presented here include Ha emitters, but the completeness estimated from 
excluding the Ha emitters is found to be the same. The separating line rep- 
resents the exact redshift range of the nan'ow-band filter when assumed to 
be a perfect top-hat function. Results are also presented in Figure [aT] 



Sample 
COSMOS 



Sources 

(Number) 
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C (%) 
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Co (%) 
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Figure Al. The completeness fraction and the fraction of sources within 
the NB filter redshift range (1— contamination) as a function of photometric 
redshift width for the COSMOS field. The dashed line indicates the narrow- 
band filter cut-off assuming a perfect top-hat filter. 



The results for different photo-z cuts can be found in Table 
|Al| and Figure [AT| These clearly reveal that adopting a photo-2 cut 
matching the full spectroscopic redshift distribution of the narrow- 
band filter (0.825 < Zph < 0.865) recovers 72 per cent of the 
total number of galaxies at those spectroscopic redshifts, missing 
a significant fraction of the population. On the other hand, while 
using a wider photo-z cut increases the completeness, it also raises 
the contamination quite considerably. 

It should be noted that in the area covered by the NBj imag- 



ing there are 22 sources spectroscopically confirmed to be within 
the naiTow-band redshift range which are missed by all the stud- 
ied photo-z cuts (11 of these are Ha emitters). These present a 
median photometric redshift of 0.73 (although some of those have 
photo-2S of 1.3 — 1.4), and are much bluer (U-K~ 0.4 — 3) than 
sources recovered by the photo-zs (U-K« 1.5 — 7.5); this indicates 
that the photometric redshift selection is more likely to miss bluer, 
less massive galaxies than redder, more massive/luminous galax- 
ies. This is mostly a consequence of the fact that blue, star-forming 
galaxies present weak 4000 A break features, combined with a sig- 
nificant [Oil] 3727 emission and Balmer breaks which, combined, 
can mimic a 4000 A at lower rest-frame wavelengths and result 
in lower photometric redshifts. In order to investigate this further, 
the completeness and contamination are computed independently 
in different stellar mass bins. The results, presented in Table |A2] 
show that the completeness increases as a function of stellar mass; 
this increase can be well approximated by a 10 per cent increase 
in completeness per 1 dex in stellar mass for the current sample - 
between 9.5 and 1 1.5 in log M©. It is also found that the contami- 
nation decreases with increasing stellar mass, by approximately 10 
per cent per dex in stellar mass. These corrections are incorporated 
when computing the star-forming fraction as a function of mass. 

For the UDS field, the lack of spectroscopic redshifts does not 
allow a detailed study of the completeness and contamination of the 
underlying samples; the completeness is therefore estimated using 
the HiZELS Hq sample and the recovery fraction of this sample as 
a function of the photometric redshift cut. In the COSMOS field, 
the completeness and contamination estimated using the Hq sam- 
ple are found to be broadly comparable to those from using the 
population as a whole. Results are presented in Table |A3l 

Completeness and contamination may vary as a function of 
environment as well, particularly because galaxies are spread rather 
uniformly in the field, but groups are clustered at certain redshifts. 
Nevertheless, currently the limited spectroscopic samples do not al- 
low to reliably estimate how completeness and contamination may 
vary with environment, although any such variation is expected to 
be weak and to largely even out in the statistics, so any possible 
trends are ignored in this paper. 



A2 The effect of varying photo-z cuts 

The choice of the underlying population for estimating local den- 
sities can in principle play a significant role in such estimation in 
a number of ways: sources at significant different redshifts (photo- 
z outliers) will introduce extra scatter in the measurements; wider 
photo-2 selected samples might pick up large structures just outside 
the narrow-band filter range and bias certain regions; more impor- 
tantly, a higher or lower number of sources could change the nor- 
malisation of the star-forming fraction. In principal, the latter issue 
can be avoided by using the completeness and contamination statis- 
tics for the underlying samples, calculated in Section Al to correct 
the projected local densities (E) and star-forming fractions (SFF) 
using: 



Ec = E X 



Co 



C 



SFFc = SFF X 



C 



Co 



(Al) 



In order to assess any potential biases or errors, corrected lo- 
cal densities have been computed using all of the previously de- 
fined underlying populations. It is found that once the corrections 
are applied, the local density estimates using different underlying 
samples correlate very well with a slope of unity. Furthermore, as 
Figure[A2|shows, the trends presented in this paper are independent 
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Table A2. A completeness and contamination study as a function of stellar mass of the photometric redshift selected samples obtained in the COSMOS field. 
The results reveal that that while completeness increases as a function of stellar mass, the contamination decreases as a function of mass. 





COSMOS 109'^ Mq lO^f -0 Mq 
Sample C (%) C (%) 


IO^O-SMq 10"-" M0 IQII-SMq 
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Figure A2. The relatively small increase of the star-forming fraction within 
the field regime and fall for denser regions is recovered independently of the 
photometric redshift cut, once corrections for coinpleteness and contami- 
nation are made. Both the underlying population and the estimated local 
densities are calculated and corrected appropriately for each photo-z sam- 
ple, but the resultant behaviour of the star-forming fractions as a function 
of local density for each determination is essentially the same. 



Table A3. A completeness study of the photometric redshift selected sam- 
ples obtained in the UDS field based on the Ha sample. Completeness is 
estimated as the fraction of Ha emitters recovered as compared with the 
total number 



Figure A3. The mass dependence of the star-forming fraction presented in 
this paper is clearly recovered independently of the choice of photometric 
redshift cut. Furthermore, correcting the star-forming fraction based on the 
mass-dependent completeness and contamination also shows that the nor- 
malisation, along with the trend, agrees very well between samples. 



bias in either or both the selection of the underlying population 
and/or the mass estimates. As Figure |A3| fully details, the mass- 
downsizing trends are robust against the choice of the underlying 
population. They are also found to be equally robust against varia- 
tions in the mass estimates (cf. Appendix B). 



Sample 
UDS 



Sources 

(Number) 



Completeness 

C(%) 



0.780 < Zph < 0.910 
0.790 < Zph < 0.900 
0.800 < Zph < 0.890 
0.810 < Zph < 0.880 
0.820 < Zph < 0.870 
0.830 < Zph < 0.860 



3502 


57 


3078 


54 


2494 


49 


1874 


42 


1502 


41 


1125 


34 



of the choice of underlying population (within the photo-z range 
used). 

Mass-downsizing trends have also been presented in this paper 
- these reveal that the fraction of star- forming galaxies decreases as 
a function of mass. It is fundamental to confirm that this is a con- 
sequence of galaxy formation and evolution, rather than a possible 



APPENDIX B: STELLAR MASS ESTIMATES 

Stellar masses were computed for all Ha emitters and galaxies 
within the underlying population as fully detailed in Section [23] 
Here, further attention is given to the uncertainties and possible 
systematic errors. 

Figure [BT] shows a comparison between the mass estimates for 
galaxies in the COSMOS field obtained in this paper and those from 
[Mobasher et al.| J2007) . There is a systematic offset from zero, but 
this is fully explained by the different Initial Mass Functions (IMF) 
used: using a Salpeter IMF, as in Mobasher et al.| ( |2007 ) results in 
assuming a mass-to-light ratio higher by a factor of 1.8 when com- 
pared to the Chabrier (2003) IMF used here; this results in a differ- 
ence of « 0.17 dex in stellar mass. The additional scatter is mostly 
due to the simpler approach used by|Mobasher et al.](|2007[l to esti- 
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Figure Bl. Comparison between the determined stellar masses and those 
from |Mobasher et 3111^2007^ and Cirasuolo et al. (in prep.). The ones com- 
puted in this study, by using a much larger number of bands and by taking a 
full SED fitting approach fixing the redshift are much more reliable than 
those estimated from rest-frame V luminosities combined with artificial 
rest-frame B-V colours by |Mobasher et al^pOO?^ (those authors also use a 
Salpeter IMF; this results in estimating stellar masses which are +0.17 dex 
higher). Nevertheless, there is no significant offset between both estimates 
(once one corrects for the different IMFs), and the scatter is fully consis- 
tent with the errors expected from the Mobasher et al determinations. As 
for the comparison with the preliminary results from Cirasuolo et al., the 
differences (scatter) mostly arise from the use of different models, different 
aperture corrections and, to a lower extend, the use of z = 0.845 in this 
paper instead of the best photo-z. 



mate stellar masses from visible luminosities and artificial colours, 
together with their poorer photometric redshifts. Figure |B1| also 
compares the mass estimates for the UDS field with those derived 
by M. Cirasuolo et al. (in preparation), revealing a lower degree 
of scatter, which mostly arises from a different aperture correction 
(M. Cirasuolo et al. use a constant aperture correction assuming all 
sources are point sources, while in this paper individual corrections 
are used; the constant aperture correction matches very well to the 
average of the individual corrections used in this paper), and dif- 
ferent stellar synthesis models (with different parameters). Overall, 
the scatter between the different determinations is relatively small 
and there are no systematic trends. 
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